Microwave controlled phase coherent transport in mesoscopic S-N-S structures 
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We show that transport through a superconducting quantum point contact biased at subgap voltages 
is strongly affected by a microwave field. The subgap current is increased by several orders of 
magnitude. Quantum interference among resonant scattering events involving photon absorption 
is reflected as an oscillating structure in the I-V curve. We also discuss how the same interference 
effect can be applied for detecting weak electromagnetic signals up to the gap frequency, and how 
it is affected by dephasing and relaxation. 



I. INTRODUCTION 

The electrical properties of small sized conductors at 
low temperatures are significantly affected by the finite 
spacing of electron energy levels, resulting in a number 
of mesoscopic phenomena. A few examples of equilib- 
rium or near equilibrium mesoscopic effects are magnetic 
moment fluctuations in quantum dots and rings as well 
as universal conductance-fluctuations in quantum point 
contacts and nanowires.cn 2 ] The opposite limit, of strong 
non-equilibrium behavior, appears for example when ex- 
ternal fields produce significant redistribution of elec- 
trons between quantized levels or modify the quantum 
states. In this limit coherent time dependent dynamics 
in mesoscopic conductors can be studied. Two possible 
approaches are nonlinear transport phenomena in strong 
electrical fields and optical interlevel transitions under 
infrared irradiation. 

An interesting possibility for coherent time dependent 
dynamics appears in mesoscopic structures where the po- 
sition of quantized electronic levels slowly (adiabatically) 
change in space or time. Under this condition external 
alternating fields will produce a strong resonant inter- 
mode coupling which will take place only in the vicinity 
of certain points in space or at certain times when the res- 
onance condition (the inter-level distance coincides with 
the energy, Ku, of the external field) is satisfiedtra. The 
adiabatic dynamics of the energy levels together with 
the well defined inter-mode coupling events results in 
a novel type of non-equilibrium mesoscopic phenomena. 
In Ref. H a ballistic channel with varying width is stud- 
ied. It is shown that interlevel transitions induced by 
infrared radiation may have a similar effect as impurity 
scattering and significantly influence transport through 
the channel. In this system, impurity scattering phe- 
nomenon like quantum localization of conducting elec- 
trons, and phenomenon foreign to impure systems like 
resonant backscattering of electrons!] will be present. As 
an example of a non-equilibrium mesoscopic phenomenon 
induced by inter-mode coupling events localized in time 
one can consider how energy is stored in an disordered 
mesoscopic ring subject to a time-dependent magnetic 



flux. The interference pattern which appears in this 
system as a result of the different temporally localized 
inter-mode transitions determines the various pictures of 
coherent energy accumulation that has been discussed 
in Ref. fj], ||. Another very interesting subject to study 
is coherent dynamics of temporally localized inter-mode 
transitions in ballistic Josephson junctions subject to in- 
frared irradiation. This is the subject to be addressed 
here. 




FIG. f . Single mode adiabatic superconducting microcon- 
striction with normal region of length L, width W and trans- 
mission coefficient D. The phase difference over the junction 
is defined as <j> = (f>2 — 4>i ■ 

It is well known that the current through a voltage- 
biased superconducting quantum point contact (SQPC) 
is carried by localized states. These states, called An- 
dreev states, are confined to the normal region of the 
contact. The energy of the states — the Andreev levels 
- exist in pairs, (one above and one under the Fermi 
level), and lie within the energy gap of the supercon- 
ductor, with positions which depend on the difference 4> 
in the phase of the superconductors across the junction, 
see Fig. [l| The applied bias affects this phase difference 
through the Josephson relation, (f> = 2eV/%. With a con- 
stant applied bias V much smaller than the gap energy 
A, cj> will increase linearly in time, and the Andreev lev- 
els will move adiabatically within the gap. This motion 
is a periodic oscillation in <j>, indicating that no energy is 
transfered to the SQPC and a pure ac current will flow 
through the contact. This is actually the ac Josephson 
effect. 

We wish to study this system in a non-equilibrium sit- 
uation; one way to accomplish this is by introducing mi- 
crowave radiation with a frequency u> ~ A/Ti, which will 
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couple the Andreev levels to each other. The radiation 
will represent a non-adiabatic perturbation of the SQPC 
system. However, if the amplitude of the electromagnetic 
field is sufficiently small, the field will not affect the adi- 
abatic dynamics of the system much unless the condi- 
tion for resonant optical intcrlcvcl transitions is fulfilled. 
Such resonances will only occur at certain moments de- 
termined by the time evolution of the Andreev level spac- 
ing. The resonances will provide a mechanism for energy 
transfer to the system to be nonzero when averaged over 
time and hence for a finite dc current through the junc- 
tion. The rate of energy transfer is in an essential way 
determinedJjftf the interference between different scatter- 
ing eventstM'cl, and therefore oscillatory features appear 
in the (dc) current-voltage characteristics of an irradi- 
ated SQPC. Dephasing and relaxation will affect the in- 
terference pattern and may even conspire to produce a 
dc current flowing in the reverse direction with respect 
to the applied voltage bias. The mechanism behind this 
negative resistance is very similar to the one responsible 
for the "somersault effect" discussed by Gorelik et aljj. 

II. THEORY 

The system we are considering is a one-mode super- 
conducting quantum point contact (SQPC) characterized 
by an arbitrary energy independent transmission coeffi- 
cient D. A_aiulti-mode junction where a single mode 
is dominantlij can also be considered. The junction is 
weakly biased, eV -C A, and a microwave field is applied 
to the normal region. The length, L, of the normal re- 
gion is taken to be much shorter than the superconductor 
coherence length £q. Further, we choose to neglect the 
proximity effect, which is acceptable in a point contact 
geometry, W « £o, where the current density is low ex- 
cept at the very center of the junction, so that the order 
parameter can be assumed to be constant on cither side 
of the contact, which is placed at x — 0. As an Ansatz 
we will work with quasiclassical envelopes u±(x,t), left 
and right movers, of the two-component wave function 
iff(x,t) — u + (x, t)e lkFX + u-(x, t)e~ tkFX , as solutions 
to theptime-dependent Bogoliubov-de Gennes equation 

(BdG)la, 

ihdu/dt = [H + V g (x, t)]u , (1) 

where u = [it + ,w_] is a four-component vector and the 
gate potential V g (x, i) = V UJ (x)a z costot induced by the 
microwave field oscillates rapidly in time. In this text we 
will use (111,112) to denote the scalar product between 
the four-component vectors. The Hamiltonian, Ho for 
the electrons in the electrodes of the point contact, 

H = — ihvF<J z T z d/dx + (2) 
A[cos(<j>(t) /2)a x + sin(<f>(t)/2)sgnx(Ty], 

where <7j and tj denote Pauli matrices in electron-hole 
space and in ± space respectively. 



The function u, which is smooth on the scale of the 
Fermi wavelength, has a discontinuity at the contact 
which is determined by the transfer matrix of the SQPC 
in the normal statp-,and is described by the following 
boundary conditioned: 

u(+0) = (1/Vd)(1-V1t 1 ,)u(-0), R=l-D, (3) 

where D is the junction transparency, < D < 1. This 
form, derived in the appendix, is adequate for our model 
and is accurate to order L/£ . 

It is well known that there is a discrete set of energy 
levels within the energy gap of the superconductors in a 
SQPC. Because of spatial quantization in the transverse 
direction, the energy spectrum of the Andreev bound 
states consists of a discrete-set of pairs of levels labeled 
by the quantum number no, 

E n ,± = ±Ay/l -D n (k F ) sin 2 (0/2). (4) 

Energy is measured from the Fermi level and the trans- 
mission coefficient D n is related to the propagation of 
normal Fermi level electrons through the microconstric- 
tion. The energy of these levels is governed by <j> — <fo — </>i 
which is the superconductor phase difference over the 
junction. 

We are considering a single- mode junction, so n = 1 
and we have two Andreev states, one above and one under 
the Fermi level, which we label as uf(x,t) and u7(x,t). 
By assuming that the voltage drop occurs only at the 
point contact we can use the Josephson relation for the 
phase difference, <j) = 2eV/h. This assumption together 
with the restriction that eV <C A, allows us to treat 
u^(x,t) and u7(x,t) as states moving adiabatically in 
time, through <f>(t) = 2eVt/%. This coupling between the 
phase difference, <f> and time, t, allows us to label the 
Andreev states as follows, u^ t ^(x) and u^, t -.(x). In this 
picture we have a pair of energy levels which will oscillate 
periodically in time with a period of T p = irh/eV, see 
Eq. ^. These states carry a current when populated and 
for a junction that is short (L <C £o)> all supercurxent 
through the constriction is carried by these stateal3. A 
study of the boundary conditions for the population is 
clearly motivated. 

A. Boundary condition at </> = 2im 

The interesting region is at <f> = 2im (see Fig. || point 
C), where the levels approach the continuum and the evo- 
lution of states cannot be considered as adiabatic, even 
if we assume a small applied voltage or a weak electro- 
magnetic field. Let us study the system in the vicinity of 
= (n — 0) and also t — 0. 

Using the Heisenberg uncertainty relation, StSE w ft 
and defining 5E = A—E(t) = ESt, then assuming Dwl, 
we approximate the duration of the non-adiabatic region 
as, 
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st = Jn/E^h(Ae 2 v 2 r 1/3 - 



(5) 



If we want St -C T p , we find the condition, eV/A <C 7r 3 
which means that we can safely treat the non-adiabatic 
region as short. 

The quantity that we are interested in is the transition 
probability between the state u + at the times t\ <C — St 
and ti ^> St. Quantum mechanically this can be ex- 
pressed as the matrix element ((uj, t U(t 2 , *i) u J( tl )))> 
where U(t2,t\) is the exact propagator with our Hamil- 
tonian, Eq. ||. 

Let us introduce the unitary symmetry operator, A = 
-Pc x r z , under which our Hamiltonian is invariant (P is 
the parity operator in x-space). This property is there- 
fore valid on both sides of a; = 0. An important point 
now is to show that the boundary condition at x = 
does not destroy this symmetry. Let us insert Au into 
the boundary condition, 

Au(+0) = (l/VD)(l - VRTy) Au(-O) (6) 
& Pa x r z u(+0) = (1/VD){1 - VRr y )P<J x T z u(-0) (7) 



apply a x T z from the left and P, 

u(-0) = (1/VD)(1 + VRt v )u(+0) 
apply from the left, (l/\/Z>)(l - VRt v ) 

{1/Vd){1 - VRTy)u(-O) = u(+0) 



(8) 



(9) 



We have recovered the boundary condition (g) which 
proves that it is invariant under A. 

Now, we can conclude that both the Hamiltonian, Ho, 
and the boundary condition at x — are invariant under 
A. This implies that at any time any non-degenerate 
eigenstate of the Hamiltonian is an eigenstate of A with 
the eigenvalue +1 or —1 and that this property persists 
during the time evolution of the state. Specifically, if we 
take a state on each side of 6 = 0, 



Au+ 



Aiu+ 
A 2 u+ 



(10) 
(11) 



Next we insert = YuJ, where the operator Y = 
a x T y changes the sign of the phase <f> into Eq. 
apply a x T y from the left and arrive at 
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and 



Auj 



A 2 uj 



(12) 



which shows that Ai = — A 2 . This means that the two 
states, and are orthogona . This is consistent 
with the results of Shumeiko et.al.H, who have shown 
that the Andreev state wave functions are 47r-periodic 
whereas the energy levels and the current are 27r-periodic. 
Since the state evolving from the adiabatic state u J( tl ) 

is orthogonal to the adiabatic state u J( t2 p the proba- 
bility for an adiabatic Andreev state to be "scattered" 



into a localized state after passing the non-adiabatic re- 
gion is identically zero. In reality, the Andreev state as 
it approaches the continuum band edge decays into the 
states of the continuum. Such a decay corresponds to 
a derealization in real space and is. the mechanism for 
transferring energy to the reservoirEj. 

The orthogonality property shown above guarantees 
that the coherent evolution of our system persists during 
only one Josephson oscillation and that the equilibrium 
population of the Andreev levels is reset at each point 
6 = 2tt?tO. This imposes the boundary condition 



6+(27rn + 0) = 0, + 0) = 1 

at the beginning of each period. 

E(4» c y y 



(13) 
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FIG. 2. Time evolution of Andreev levels (full lines) in the 
energy gap of a gated voltage-biased, single-mode SQPC (see 
inset). A weak microwave field induces resonant transitions 
(wavy lines) between the levels at points A and B and the 
level above the Fermi energy becomes partly occupied to an 
extent determined by interference between the two transition 
amplitudes. Non-adiabatic interactions release the energy of 
quasiparticles in the (partly) occupied Andreev level into the 
continuum at point C, where the Andreev states and the con- 
tinuum merge into each other (represented by dashed arrows, 
see text) and the initial conditions for the Andreev level pop- 
ulations are reset (filled and empty circles). 



B. Effect of radiation on the population 

Now we know that we have a system which is periodic, 
and it is possible to study the effect of the applied high 
frequency radiation on the population of the Andreev 
states and the dc-current. We will work with applied high 
frequency radiation that will induce transitions between 
the two Andreev levels, described by the eigenfunctions 
to the Hamiltonian, u + and u~ . The high frequency field 
will radically effect the system when the condition for 
resonance is fulfilled. This means that even at small am- 
plitudes, V w <§; A, when the energy of the field matches 
the interlevel distance between the two Andreev states, 
E + (t) — E-(t) = Huj or 2E + (t) = Hm, there will be a large 
change in their population. In this situation we cannot 
use perturbation theory and instead we will use the res- 
onance approximation, where u is taken to be a mixture 
of the two states u + and u~ , 



u(x, t) 



>(*) 



(x)e 



-i(uj/2)t 



(x)e 



(14) 
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Inserting this Ansatz into the BdG equation (]l]) and aver- 
aging over fast oscillations we arrive at the following cou- 
pled differential equation for the population coefficients, 
b=[b+,b-}; 



ib(t) 



vl_/n -Su+,-(t) 



b(t), 



(15) 



where 5w+ t -(t) = [2E + (t) — huj]/2h, is a mea- 
sure of the deviation from resonance. Vi\ = 

dxV UJ (x)/2(u + ,(7 z u _ ) is the matrix element between 
the Andreev states.Q 

The matrix above is very close to model Hamiltonians 
uscdrta-study Landau- Zener dynamics of two level sys- 
temsta'ta. Before we start to study the dynamics of this 
equation, we need to find an expression for the current 
through the SQPC. 



III. CURRENT 

The current carried by populated Andreev levels in 
equilibrium can be calculated from / = 2e/hdE n ((f>)/d<fi. 
This equation is also valid for adiabatically moving An- 
dreev levels, however, it needs modification when inter- 
action with a high frequency field is considered. 

We start with the equation for continuity, where p e — 
e(u, a z u) is the charge density in a normal conductor 
(A = 0): 



d t Pe + d x I = 0. 



(16) 



This provides us with the equation for the current, 
I{x,t) = —evF{u,T z u). In superconductors, the charge 
conservation equation takes form: 



dtPe + d x I = 1 
2eA 

I = — — [sin (</>/2)sgn(x)(u, a x u)- 
cos (0/2) (u, <Jy\i)\ . 



(17) 
(18) 



By multiplying this equation with sgn(x) and integrat- 
ing over the whole x-axis we obtain the current at the 
junction (x = 0) as 



/oo 
dxsgn(x)(l-d t p). 
-OO 



(19) 



The current Idc averaged over one period of Josephson 
oscillations reads 



l dc 



dx sgri(x) 



Jdt 



P\o P 



(20) 



The first term in this equation characterizes the charge 
exchange between Andreev states and condensate. The 
second term describes the ac-component of the current 
which is connected with a periodic charge redistribution 



in the vicinity of contact. This term is equal to zero 
which is shown in Appendix B. After a transformation of 
the quantity I (see Appendix B), Eq. |2(] can be presented 
in the equivalent form, 



I, 



dc 



2e 



dt 



dx 



iti 



d(u, u) 
dt 



V g (u, a z u) 



(21) 



To find an expression which is valid for our situation with 
the applied field and with only two energy levels we can 
insert the resonance approximation, Eq. \L% 



2e , . 
I dc = — A- 
Tvk 



\b+(t = T p )f 



(22) 



A detailed derivation can be found in Appendix. 

To check the validity of the current expression, let us 
study energy conservation. We have two sources of en- 
ergy, the applied field and the applied bias. If we con- 
sider a single Josephson period, the energy absorbed by 
the system from the voltage source is I ■ V ■ T p and from 
the applied field it is hu>\b + (T p )\ 2 . Energy conservation 



for the period can be stated as, E Q 



E„ 



rbed ■ 



the energy leaving the system stated as 2A|6 + (T P ) 

2A|6+(T P )| 2 = I ■ V ■ T p + hu;\b + (T p )\ 2 , 
the current can easily be found as, 



I = 



1 



■{2A\b + (T p )\ 2 -hLj\b + (T p ) 



or, inserting the value of T p , 



i<A-*V(T,)|. 



with 



(23) 



(24) 



(25) 



This is exactly the same result as the more tedious calcu- 
lation found in the appendix. So, the current can be un- 
derstood as voltage bias mediated, photon assisted tun- 
neling. 

To sum up, we have shown that the dc-current through 
a SQPC is linearly dependent on the population of the 
upper Andreev state at the end of the Josephson oscil- 
lation period. This follows from the periodicity of the 
population of the Andreev states, Eq. |l3j. In equilibrium 
without an external field present to shift the population 
of the Andreev states, which, at T = 0, is only in the 
lower state, there will be no current. The introduction 
of an external field, which is considered here, affects the 
population of the Andreev states in a non-trivial man- 
ner and therefore also the current. The dynamics of the 
Andreev state populations under radiation, governed by 
Eq. p~3, is therefore the next subject. 



IV. DYNAMICS OF THE ANDREEV STATES 
UNDER RADIATION 

The current through the SQPC depends linearly on 
the population of the upper Andreev level at 4> = 2nn~ , 
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the end of the Josephson period. Depending on the fre- 
quency of the applied field there will be either one, two 
overlapping, or two well separated resonant events dur- 
ing one period. The resonances will be Landau-Zener like 
since the energy levels move in and out of resonance peri- 
odically in time. For a systematic study of the dynamics 
we will first look at (1) the case of a single resonance and 
overlapping resonances, then we will consider (2) the case 
of well separated resonances. 

Initially we can discuss two different limits to qualita- 
tively understand the dynamics. The first is when the 
applied bias is large. In this limit the energy levels oscil- 
late at a high rate and there is no time for a resonance to 
shift the population, the population will be untouched. 
The opposite limit is when the applied bias is very small, 
then the energy levels will oscillate very slowly. In this 
limit the upper level will be fully populated and then 
fully depopulated during one period, once again the pop- 
ulation will be untouched. When the population is un- 
touched there will be no dc current, and we can conclude 
that in these limits the dc current will be small, and that 
the interesting region will be for intermediate bias. The 
same discussion can be applied for the limits of a weak 
and a strong applied external field. 

A. Single and overlapping resonances 

Consider the case when the applied field matches the 
interlevel distance in the middle of the period, at <f> = n, 
see inset of Fig. |^a. At this frequency there is only one 
resonant event. We wish to sweep over a large range 
of applied bias in order to calculate the current voltage 
characteristics. The term that governs the dynamics is 
the distance from resonance, 

5u)+-(t) = |(2Aa/i- DJk F )sin 2 (cf)l2) - hto), (26) 
n v 

which has a parabolic time dependence for <f> sa tt. This 
means that we need to solve the system numerically by 
time iteration of the differential equations (|l5|). 

Numerically it is also possible to study the behavior 
of the current when we first have a single resonance and 
then in small steps increase the frequency of the applied 
field producing two resonances, Fig. |^. The onset of the 
interference between the two possible paths to the upper 
levels can clearly be seen. The origin of this interference 
is that the phase collected by a quasiparticle depends on 
which level it populates between the two resonant events. 
This can be compared to the classic double slit experi- 
ment in optics. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 

eV/A »i 



0.35 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 

eV/A x, 



FIG. 3. Numerically calculated, radiation induced current 
I c = eAD/h. Top, a single resonance and bottom, small de- 
viation from single resonance (oj/(tj S i ng i e = 1.0155) and the 
appearance of interference. The time (position) of resonance 
is shown in the inset figure. In this case the differential equa- 
tions describing the system have been solved numerically since 
the time dependence of the energy levels is parabolic and the 
two resonances can effect each other. In both plots A = 2K 
and V+-/A ta 0.0036. 

When the strength of the applied field is increased, 
new oscillations are added to the current, and the enve- 
lope of the current is shifted to higher bias, see the right 
part of Fig. ||. We interpret these "new" oscillations as 
a precursor to Rabi oscillations. At <fi ~ the time de- 
pendence is close to zero, dE/dt ~ 0, and an analogy can 
be drawn to a "long channel" where Rabi oscillations ap- 
pear between two energy levels with-a. period depending 
on the strength of the applied fielded. The population 
is pumped up and down a number of times during the 
resonance period resulting in the extra oscillations found 
in Fig. | 
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eV/A 



I, /1 0.15 

dc c 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 

eV/A 

FIG. 4. Plot of the single resonance current for an in- 
creasing strength of the applied field. In the upper plot, 
V+.-/A « 0.036 and the lower V+.-/A ft! 0.12. In both 
plots A = 3K. 



B. Separate resonances 

For the case when there are two separate resonances a 
linear approximation can be applied to the time evolu- 
tion of the energy leizsls included in Eq. This is the 
Landau-Zener modelEU and it is valid with the restriction 
that the applied field has to be weak. Then an asymp- 
totic solution to the differential equations is available, 



1 



and 



1 



7 = ir\V+^\ 2 /\dE/dt\, 



(27) 
(28) 



where d 2 is the probability of changing state after one 
resonance, and r 2 is the probability of being "reflected" 
back to the same level. The first step is to construct a 
unitary scattering matrix for each resonance, Sa and Sb, 
using symmetry arguments which are valid for our sys- 
tem, Eq.([l5|), like & = a z S a z and Sa,b = OxSb,a<Jx 
we have 



S A = 



r 
de 10 



-de 
r 



-iB 



, Sr — 



r 

-de~ iB 



de ie 
r 



(29) 



where 9 is the phase that is transferred by a transi- 
tion. To take care of the phase gathered between the 
transitions ($), we also introduce a propagation matrix 
U = cxp(— iOzQ) with 



$ = 



1 



2eV 



(30) 



The phase $ is found through the coupled differential 
Eqns. (plf) with = 0, the solution far away from 



resonance, 4>a and <Pb being the two events of resonance, 
see Fig. |[ 

By applying the scattering matrices, Eq. (|2S|), and the 
propagation matrix U to the initial state bo — (0, 1) in 
this fashion b(T p ) — Sb U Sa^Oj we arrive at 



\b + (t = T p )\ 2 = 4r 2 d 2 sin 2 (9-$). 



(31) 



I. II 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 



eV/A 

FIG. 5. Calculated dc-current, (I c — eAD/h) when the 
frequency of the applied radiation is such that there are two 
separate, non-overlapping resonances between the Andreev 
levels, see Fig. 2. In this case the resonances are placed such 
that the time variation of the energy levels can be approxi- 
mated as linear and that the two resonances can be treated 
as separate. V+-/A 0.048 and A = 3K. 

The current appears in the subgap region with the fol- 
lowing form, and using Eq. we acquire, 



8e 



d 2 (A-^)sin 2 (9-$). (32) 



The sine term in Eq. (|3^) represents the interference of 
the two probability amplitudes representing the different 
paths to the upper level. The phase collected by a parti- 
cle will be different depending on which level it populates 
between the two resonances. When the distance between 
the two resonances is increased the number of oscillations 
will also increase, see Fig. || and compare with Fig. ||b. 
This onset of this interference effect can be studied by 
comparing Figures ^ and |^. The origin is simple and 
presents a new method of investigating the nature of An- 
dreev states. The details of this interference is studied 
and presented in Ref. |[ 
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FIG. 6. Plot of the current as a function of inverse volt- 
age, for a well separated resonance, Tiuj /2A ~ 0.76, with an 
increasing strength of the applied field. In the upper plot, 
V+.-/A « 0.024 and the lower V+.-/A « 0.12. In both plots 
A = 3K. 

When the strength of the applied field is increased such 

that V-\ /A oc 0.1, and \b + \ 2 is calculated numerically 

through Eq. [l5[ we find additional oscillations, which the 
analytical approach failed to produce. They modulate 
the 1/V periodic oscillations, which have their origin in 
the systems characteristic period T p = hir/eV. See the 
left part of Fig. |^ where the inverse current is plotted for 
three different values of V-\ 

When an effort to study the mechanism presented here 
experimentally is undertaken, the most probable method 
should be to slowly increase the strength of the applied 
field and study the FV curve. This can be done for sev- 
eral different regimes, (1) for a low or high applied bias 
and (2) for different frequencies of the applied field, rang- 
ing from a single resonance to two well separated ones. 

In Fig. ^ two plots are presented for different values 
of the applied bias when we have only one single reso- 
nance. For the case when eV/A = 0.012 there appears 
weak oscillations for higher V±. These can probably be 
understood as Rabi-like oscillations, we can approximate 
the interlevel distance as constant during the resonance 
and the population can be pumped up and down several 
times. 
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+,- +.- 

FIG. 7. The current as a function of the applied field 
strength for the case of a single resonance. In the left fig- 
ure eV/A » 0.012 and in the right figure, eV/A ft! 0.12 and 
in both A = 3K. 

Another case that can be studied is the current as a 
function of the applied field's frequency, tku. The familiar 
oscillation will appear even in this case, see Fig. ||. By 
slowly increasing the frequency of the applied field and 
thus shifting the position of the resonances, the sine term 
in Eq. B3 is modulated. 
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co/2A co/2A 
FIG. 8. The current as a function of the frequency of 
the applied field. In the left figure eV/A « 0.019 and 
V+,_/A w 0.036, while in the right figure, eV/A « 0.12 and 
V+-/A « 0.12 In both plots A = 3K. 



C. Relaxation and dephasing 



Up to this point we have studied an ideal system, as- 
suming that Landau-Zener like transitions dominate. We 
have neglected the effect of voltage fluctuations and in- 
teractions with external degrees of freedom. If we as- 
sume that the relaxation and dephasing times, r, are 
short compared to the duration of the non-adiabatic res- 
onances, St, Eq. |, we can phenomenologically study the 



7 



effect of such interactions on our system. The techniquecHI 
we will use models dissipation by adding a term to the 
time evolution equations of the density matrix for our 
two-level system. The equations have the following form, 



Pnn(t)=--[H (t),p(t)} n 

Mt) = -~[H (t),p(t)], 



Pnn Pn n 



Pnn' 



, n^n', 



(33) 
(34) 



with r as the characteristic time for relaxation of the 
system to the equilibrium population, p eq , and as the 
characteristic time for dephasing. 

The exact form of the density matrix for the system 
considered here will be a 2 x 2 matrix for the discrete two 
level system of the Andreev states (To avoid confusion we 
will use cti, (72, er 3 to denote the Pauli matrices in this 
discrete space). The Hamiltonian for the density matrix 
is found in Eq. [l5|., and the diagonal elements p\\ and 
P22 will represent the population if the upper and lower 
Andreev levels. The resonant events which are unaffected 
by the dissipation can still be modeled with the scattering 
matrices introduced earlier, Eq. E9, 



P 



Sa,b P S a 



A,B ' 



(35) 



To demonstrate the effect of relaxation and dephasing 
we will calculate the dc-current for well separated reso- 
nances. In this case the system is strongly effected by 
the high frequency field only in the vicinity of t = ti : 2- 
Therefore during almost the whole interval of time, T p , 
one can neglect V g in Eq.(l) and calculate the instanta- 
neous current I(t) by using the expression 



(36) 



This expression is an extension of Eq.(20) into the case 
of relaxation. As a result in adiabatic limit h<p « A, 
the dc-current is given by 



T p n 

1 



Tr(a 3 p)dE/dcj) dt 



9E+(t) 
dt 



[2p u - 1] dt, 



(37) 



2eV/h 



where we have used the Josephson relation < 
and that Tr{p) = 1. 

By solving the time evolution of the density matrix for 
the whole Josephson period and imposing the boundary 
condition p eq = (1 — o- 3 )/2 at t = 0, we arrive at the 
following, 



Pu = 





d 2 e -(t-tx)/r 
f(r 2 ,r^r) e-(*- 



■t 2 )/r 



, o < t < h 

,t 1 <t<t 2 

,h<t<T p 



(38) 



f(r 2 , T4> ,T) = \d 2 (l - e -(**-ti)/T) + 2r 2 d 2 (e-( t ^ tl )/ T - 



-(ta-tl)/T> 



cos2(6 + $)) 



(39) 



where d 2 = 1 
we find, 

Idc(V) 



Inserting these expressions into Eq. |37] 



VT V 



d 2 huj 



-<.u- tl )/T_A + 



d 2 r 2 

— / E+(t)dt+ 
T Jti 



f(r 2 ,r,,r)^Ae-^-^ 



E + (t)e-^- t2 ^ T dt 



(40) 



This formula falls back into the previous current ex- 
pression that we have derived, Eq.^2| when — + oo and 
r — > oo. From expression ( p9|) one can see that the os- 
cillations of the current exponentially decrease when the 
dephasing time becomes smaller than T p . At the same 
time when the relaxation time r is decreased the cur- 
rent drops tOp-aero as r. The I-V curves for the case 
T p < sa 2r E3 are presented in Fig. ||. 

A curious effect occurs when the relaxation time is such 
that the upper level is fully depopulated when tf> = ir, 
in the middle of the Josephson period. The result is a 
negative current when the dc current is calculated for 
a positive bias, our system could under special circum- 
stances show a negative conductance see Fig. ^. The rea- 
son is that the upper level is mainly populated when it's 
derivative is negative and then the lower level is highly 
populated for the part of the period when it's derivative 
is negative. The effect should be most pronounced when 
the energy of the applied field, Hco, approaches the energy 
of the gap, 2A. 

Actually, the physical mechanism which is behind the 
appearance of this negative resistance is similar one dis- 
cussed irD. 



I/I_c 




100 



with 



8 



FIG. 9. The dc current calculated with dephasing and 
dissipation phenomenologically added. In both cases 
~ r ps 4 ■ 10~ 10 s. The frequency of the applied field is 
such that 7kj/2A as 0.75 in the left plot hu/2A « 0.95 in the 
right. Note the fact that the system shows a negative conduc- 
tance in the right plot. The amplitude of the applied field is 

such that \V-{ |/A « 0.025 in both plots, and for reference 

the case when r, T0 — > oo are included (dotted curves) . 



APPENDIX A: BOUNDARY CONDITION AT 

X = 



V. CONCLUSIONS 



We have shown that irradiation of a voltage-biased 
superconducting quantum point contact at frequencies 
lo r*> A can remove the suppression of subgap dc transport 
through Andreev levels. Quantum interference among 
resonant scattering events can be used for microwave 
spectroscopy of the Andreev levels. 

The same interference effect can also be applied for 
detecting weak electromagnetic signals up to the gap 
frequency. Due to the resonant character of the phe- 
nomenon, the current response is proportional to the ra- 
tio between the amplitude of the applied field and the 
applied voltage, / ~ |V±| 2 /AeV. At the same time, for 
common SIS detectors a non-resonant current response 
is proportional to the ratio between the amplitude and 
the frequency of the applied radiatior£3), I ~ |V±/cj| 2 , 
i.e. it depends entirely on the parameters of the external 
signal and cannot be improved. 

A large part of the text is devoted to a thorough study 
of the dynamics of adiabatic Andreev levels in a biased 
SQPC subject to microwave radiation. We have pre- 
sented calculations of the dc-current where we vary either 
the frequency or the strength of the microwave field, or 
the voltage bias. These results are open for experimental 
investigation since a shift in the population of the An- 
dreev states is directly reflected in the current through 
the SQPC. 

Finally, we note that the classic double-slit interference 
experiment, where two spatially separated trajectories 
combine to form an interference pattern, clearly demon- 
strates the wave-like nature of electron propagation. For 
a 0-dimensional system, with no spatial structure, we 
have shown that a completely analogous interference phe- 
nomenon may occur between two distinct trajectories in 
the temporal evolution of a quantum system. 
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FIG. 10. Schematic for SNS junction with scatterer. 

Here we are going to calculate a boundary condition 
for u± (x, t) at the normal region with a scatterer. This 
is done under the assumption that we have one scatterer 
located at x — and that A -C fi. 

In the left superconductor we make the Ansatz, 

* 1 = ( bl ) ^ + ( d\ ) ^ (A1) 
and analogously for the right superconductor, 

^ = (Z) elkx+ {Z) e ~ lkx - (A2) 



The Ansatz for the normal region will be two linear 
combinations of the possible scattering states, see fig. [ll], 




FIG. 11. Scattering states, from left and right. 
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and for x > 0: 







N 



At + Bip \ 
Cr + Dip ) 



-ikx 



(A5) 



By matching through this structure, Fig. [LO] with the 
Ansatzcs as defined above, we arrive at the following con- 
dition, 









= l/r| 









1 -ip 

1 -i/o 

ip 1 

ip 1 





/oa\ 


) 


b 2 




c 2 







(A6) 



Using the notation with Pauli matrices we can write this 
boundary condition as, 



9(-x) = -±=(l+T v y/R)9{x) 



V(x) 



(A7) 



(A8) 



APPENDIX B: CURRENT EXPRESSION FOR A 
BIASED JUNCTION 

In this appendix we derive an expression for the time- 
dependent current through an SNS structure when a 
weak bias is applied. An applied bias forces the An- 
dreev levels to oscillate introducing a time-dependence, 
which makes the validity of the commonly used expres- 
sion, I — (2e/h)dE/d(p, questionable. 

We wish to find an expression for the current that 
includes time-dependence. The time-dependence enters 
both through the phase difference tj>(t) — 2eVt/h and the 
applied field V g (t). 

According to Eq. [l8[ charge conservation in the super- 
conductor reads: 

h{d x I + d t p e ) — 2eA [sin (</>/2)sgn(x)(u, a x u)— (Bl) 

COS (0/2)(u,(T a u)] , 

Using the identity 

ind t (n,u) = (u,Hu), (B2) 
which directly follows from the BdG equation (1), we find 



i(o,t) = |r 1 



dx 



. <9(u,u) • 

in — ^ V g (u, a z u) 



^d t (^J Pdx-J^ Pdx^j. (B3) 



To calculate the charge density p e , one needs to exam- 
ine the properties of the Andreev state wave functions. 
This is done by evaluating the expression for the charge 



density p e = e(u, a z u) with the resonance approxima- 
tion Eq. ( [l4| ) inserted, the result is a linear combina- 
tion of the terms (u , <r z u ) and (u , <r a u T ). It follows 
from Eq.(||) that the 4- vectors u T satisfy the symme- 
try equation = r z a x (tU) . This immediately gives 
that (u ± ,cr^u ± ) = — (u ± ,cr 2 u ± )* = 0. Furthermore , at 
\E\ = A (t = 0,T P ), the Andreev state wave function in 
a contact with normal electron scattering (R ^ 0) has a 
form similar to the wave functions in a perfect constric- 
tion (R = 0) (Rcf[7]). or the case of a perfect constric- 
tion, the Andreev state wave functions u + and u~, are 
eigenf unctions of t z with different eigenvalues (±1). This 
further gives that (u± ^uT ) = 



Idc 



or 



Ldc 



Therefore the appropriate current expression is, 
l(0,t)dt = 



o 

2e ■ 1 



KT, 



V 



dt 



dx 



d(u, u) 
1 dt 



- Vg(U, (7 Z U) 



(B4) 
, (B5) 



2e 



^((u,u)lo P - / dt(V g (u,a z u)) 
o 



(B6) 



where (..) symbolizes integration over the whole x-axis. 



APPENDIX C: CURRENT EXPRESSION UNDER 
IRRADIATION 

The result in the previous appendix is a general ex- 
pression for the dc current through a short biased SQPC. 
Here we will specialize the current expression to the case 
of an irradiated single mode junction. 

For the first part ((u, u))| p we are only interested in 
the end points of the cycle, t = and t = T p . At these two 
occasions the wave functions are well defined, a mixture 
of the two available states which are orthogonal. 



u = b + (t)u + (x)e~'-!r- +b-(t)u-(x)e-'-Tr- (CI) 

We also need the time derivate, u, where we can neglect 
6+ and b~ since the population is stable far from reso- 
nance. 



u = -^b + (t)u+(x)e-'-^ 1 - l -^Lb-(t)u-(x)e- 1 ^ 1 (C2) 



iE^ 



The following is known; at t — and t = T p : Ef = A 
and Ei = —A, and that at t = 0: b~ = 1 and 6+ = 0. So 
we have, 

(u,u) = -^(|&+| 2 |u+| 2 -|6-|>-| 2 ) (C3) 

and 



<(u,u)) = -^(|6+i 2 -n 2 ) 



(C4) 
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or, using the normalization, |6 + | 2 + \ b 
((u,u))=-^(2|6+| 2 -l). 

Now we can express the first part of the current as, 

ie , iA 



flUJ 



(C5) 



(C6) 



7T tl 

— (2|6+(t = T p )| 2 -l-(2|6+(i = 0)| 2 -l) = 

2 J^(\b+(t = T p )f-\b+(t = 0)f) 

The second part is a little more complicated, here we have 
to use the resonance approximation for the wave function 
and insert the applied potential V g (x, t) = Va,(x) cos (tut). 
Then we find, 



(V g (u,a z u)} = —d t {tfa 3 b) 



the second part of the current expression now follows 



(C13) 



e 

hit 

eu> 

2n 

eui 

7T 



hir/eV 



(Vg(n,a z u)) = — 



eoj 

2 



(\b+(t = T p )f-\b+(t = 0)\ 2 ) 



e rhir/eV 




tin Jo 




2 


-\b-?)\l P = 


eoj 




2> 


= o)| 2 ). 





-^d t ($a 3 b) = 



(C14) 



Putting the two parts together, C6-C14, the current 
will be, 

/*:(* = 0) = ^ (A - ^) [|6+(t = T p )| 2 - |6+(t = 0)| 2 ] 



(C7) 



(V g (x,t)(u,a z u)) = 

[Omitting oscillating terms and defining 
V + - = (^(u+,a z u-))} = 

-(6-y+_(&+) t - (b-)Wl_b+) = 
i 

to (Re [V+ -} 6 f cr 2 b + lm[V+- ] $ <j x b) , 



where b = [b + , b~]. 

Since we need to integrate this expression over one cy- 
cle, we need more information. This can be found from 
the coupled differential equations governing the behavior 
of b. If we rewrite these equations using Re[] and Im[] 
we get the following form 



Re[V+ 



U-^W. (C8) 



ib = G 3 b 

2 



Next, if we operate from the right with 6 '03 and on the 
conjugate of the above equation with ba 3 we get the fol- 
lowing two equations, 

i$ a j = S ^u 5fcWv 2 & + i^fel^w (eg) 



and 



2 h h 

Finally, take the difference between the two, 



iSW. (CIO) 12 



i(pa z b + ba 3 b ) = ^(Re[V + _pcr 2 b + Im[V r + _]6 t (7i&). (Cll 
n 



Then we have on the right hand side, part of Eq. C7 and 
combining these two we obtain, 



id t ($a 3 b) = ^L(V g (u,a z u)) 
Two 



(C12) 
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